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Abstract. We investigate the conjectured ground state eigenvector of the 8- vertex model 
inhomogeneous transfer matrix on its combinatorial line, i.e., at rj = vr/3, where it acquires 
a particularly simple form. We compute the partition function of the model on an infinite 
cylinder with certain restrictions on the inhomogeneities, and taking the homogeneous limit, 
we obtain an expression for the squared norm of the ground state of the XYZ spin chain as 
a solution of a differential recurrence relation. 



1. Introduction 

Cuf The purpose of this article is to investigate the inhomogeneous eight-vertex model on a 

particular one- dimensional family of the globally defined parameters of the model, namely, 
with the conventions of Baxter ['^j, when rj = ir/S. More precisely, we study a certain eigen- 
vector (conjecturally, the ground state eigenvector in an appropriate range of parameters) 
of the transfer matrix of this model with periodic boundary conditions and an odd number 
of sites. Ultimately, the goal is to compare with some observations and conjectures [5, 24] 
made for the homogeneous eight-vertex model and the closely related XYZ spin chain, but 
the introduction of inhomogeneities (spectral parameters) turns out to be quite useful, as 
was previously found for the six- vertex model [11, 15, 25], a special case of the eight- vertex 
model. 

In this section we briefly describe the model and some conjectured properties aX 1] = n/S. 
The rest of this paper is devoted to showing how some of these properties arise from special- 
izing formulae for the inhomogeneous model. The main object of study will be the "partition 
function" of the model on an infinite cylinder (equivalently, a quadratic functional of the 
ground state eigenvector), for which we derive an inhomogeneous sum rule (with a certain 
restriction on the inhomogeneities, which we call "half-specialization" ) and a detailed discus- 
sion of its homogeneous limit. Note that this paper is not meant to be fully mathematically 
rigorous; firstly, it is based on a conjecture (Conj. 1) which we hope to prove in future work 
[31]. Secondly, some calculations involving theta and elhptic functions are skipped; though 
they are in principle elementary, they can be quite tedious. 

It should be noted that a special case of the eight- vertex model on its combinatorial line, 
namely the six-vertex model at A = —1/2, is much better understood [28, 1, 22, 15, 25], 
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and in this case many formulae of this work are already known and proved; we provide in 
appendix A the connection to earlier work by taking the limit to the six- vertex point. 



1.1. Inhomogeneous eight- vertex transfer matrix. The eight- vertex model is a two- 
dimensional statistical lattice model defined on the square lattice by the assignment of arrows 
to each edge of the lattice, according to eight possible local configurations around a vertex: 




R{x) 



They are given Boltzmann weights denoted by a, b, c, d which are parameterized as follows: 

a{x) = M'^V,P'')M^,P'')M^ + 2?7,p') 
b(x) = M'^V,P^)M^,P'')M^ + 2?7,p') 
c{x) = A{'^ri,p'')M^,p'')M^ + 2?7,p') 
d{x) = ^i{2r),p^)A{x,p^)A{x + 2ri,p^) 

where x is the spectral parameter and p — e^'"'^, Im r > 0, is the elliptic nome. The weights 
have period 27i and pseudo-period 2'kt, i.e., they are multiplied by a common factor when x 
is replaced with x + 2'kt. 

Ordering the edge states as (t, 4-) and (^,-<— ), these weights can be encoded into the 
i?-matrix 

/a{x) d{x)\ 
b{x) c{x) 
clx) blx) 
\d{x) a{x) J 

We shall also need in what follows the i?-matrix defined as 

fa{x) 

c{x) b{x) 

blx) clx) 
\d{x) 

where V permutes factors of the tensor product. 

The Boltzmann weights satisfy the Yang-Baxter equation and unitarity equation; in terms 

of R, these are expressed as 

(2) Ri^i+i{x)Ri+i^i+2{x + y)Ri^i+i{y) = Ri+i^i+2{y)Ri,i+i{^ + y)Ri+i,i+2{x) 
and 

(3) R{x)R{-x) = r{x)r{-x) 1 



R{x) = VR{x) 



d{x)\ 



a{x)) 
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where r(x) = "(94(0; — 2rj]p'^)'di{x — 2ri;p'^). 

We consider here the model in size L with periodic boundary conditions in the horizontal 
direction, i.e., with the geometry of a cylinder of width L. The state of the L vertical edges 
at same height on the cylinder are encoded by a sequence in {t, i}'^- The transfer matrix 
is a 2^ X 2^ matrix, or equivalently an operator on (C^)®^ with its standard basis indexed 
by {t, l}^, describing the transition from one row of vertical edges to the next; the fully 
inhomogeneous transfer matrix has the formal expression 

Tl{u\xi, ...,xl) = Tro -Roi(a;i - u)Rq2{x2 - u) . ..Rql{xl - u) 

where we use the following convention: the indices of operators R (and all other local op- 
erators) are the spaces on which they act in the tensor product . u,xi, . . . ,xl are 
spectral parameters of the model. The system has rotational invariance in the sense that 
shifting cyclically sites in the tensor product and spectral parameters leaves invariant. 
In what follows, all indices in {1, . . . , L} must be understood modulo L. 

Finally we need Pauli matrices a^'y^^^ which are local operators acting on one site; we give 
alternate names to two of them. The flip operator [a^ Pauli matrix) is F = ( ^ J ) and the 
spin operator (cr^ Pauli matrix) is a = ( g -^i )• Finally cr^ = ( q ). Denote = Y[i=i ^i- 

The transfer matrix is invariant by reversal of all spins, i.e., [T^, F*] = 0. 

1.2. Combinatorial line. In all the rest of this paper, we assume that L is an odd number, 
L = 2n + 1, and that rj = 7r/3. This second condition is what we call "combinatorial line", 
because of the occurrence of integer numbers in the ground state, as we shall see below. The 
value r] = ir/S was first noticed to have special significance by Baxter [■)]; the importance of 
odd L was emphasized by Stroganov [27]. More recently, Razumov and Stroganov [21] and 
Bazhanov and Mangazeev [4, 5] studied the model with such conditions. 

Although the work [ '] is mostly concerned with the homogeneous limit (see below), the 
following conjecture is made there (translated into our present conventions): the transfer 
matrix Tl{u\xi, . . . ,xl) possesses the eigenvalue 

L 

tL{u\xi, . . .,Xl) = Y\_i^i^i -U) + h{Xi - U)) 
i=l 

In fact, this eigenvalue is found to be doubly degenerate; in [24, 4, 5], this degeneracy is lifted 
by fixing the parity of the number of t in the eigenvector. Here we find it more convenient 
to choose a different convention, which is to diagonalize simultaneously F^. 

Note the identities at 77 = vr/S: 

r{x) = 'di{0,p'^)i!}i{x + ri,p^)-di{x + ?7,p^) = a{x) + h{x) 

1.3. Homogeneous limit. If we assume that all Xi are equal (homogeneous situation), then 
the transfer matrix commutes with the XYZ Hamiltonian, which can be written as 

1 ^ 

i=l 
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The numbering of the couphng constants will be explained later. The value rj = tt/3 implies 
that up to normalization, the three coupling constants can be expressed in terms of a single 
quantity, which we choose to be 

If we choose r purely imaginary, then as p goes from to 1, C goes from to 1. 
The coupling constants are given up to overall normalization by 

J? = — — = : Jd 



2 1+C 1 -C 

The XXZ Hamiltonian (corresponding to the six- vertex transfer matrix), is the case C = 
(or p = 0). This case was already studied in detail, as mentioned in the introduction. 

Another special case is C — ^ 1 (or p — )■ 1), for which after rescaling the weights, J2 = J3 = 0, 
so the model becomes the Ising model, but with a a^a^ interaction. The ground state 
becomes of course trivial; some details are provided in appendix B. 

The simple eigenvalue of the eight- vertex transfer matrix translates into a simple eigenvalue 
of the Hamiltonian Hi, namely 

El = -^(^2 + ^3 + ^4) 

It is conjectured to be the ground state eigenvalue of Hi. 

Many remarkable observations were made on the corresonding eigenvector, in [24, 5]. 
Its entries can be chosen to be polynomials in (, and the form of some of these polynomials 
was conjectured. We shall not discuss these conjectures here. The values at ( = (XXZ 
model) of these polynomials were calculated in [25]. 

In both [2-1, 5], the squared norm of \E'l was introduced: 
(4) \^l\'= *L 

where the normalization of the components is chosen so that they are coprime polynomials 
in (, and ^^.-[^ . . . -[-^ . . . ^(C = 0) = 1. An expression for I^E'lP was given in [o] in terms of 

certain polynomials, themselves defined by differential recurrence relations which are special 
cases of certain Backlund transformations for Painleve VI. Since the formulae are rather 
complicated, we shall not write them out here and derive our own (similar) formulae by 
specializing our inhomogeneous expressions. 

2. Properties of the ground state eigenvector 

We consider once again the eigenvector equation in size L = 2r;, + 1 

(5a) Tl{u\xi, . . . , Xi)^L(xi, ...,xl)= tiiulxi, . . . , xl)^l(xi, ...,xl) 

(5b) F,^l{xi, ...,xl) = i-ir^Lixi, ...,xl) 
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for the inhomogeneous eight-vertex transfer matrix, where we recall that tL{u\xi, . . . ,Xl) = 

Y[i=i''^i^i ~ '^)-> ^{^) = '^(^) + ^(^) = '(94(0,p^)'(9i(x + ?7,p^)'i94(x + rj,p^). The choice of 
eigenvalue of F^, will turn out convenient in what follows. 

2.1. Pseudo-periodicity. Based on extensive study of the ground state entries by computer 
for small sizes L = 3, 5, 7, the following conjecture seems valid: 

Conjecture 1. The eigenvector equations (5) possess a solution \E'i(xi, . . . , x^) whose en- 
tries are theta functions of degree L — 1 = 2n and nome in each variable Xi (generically 
non zero and without common factor); i.e., they are holomorphic functions with pseudo- 
periodicity property: 

(6a) . . , a;. + 27rr, . . .) = p-^";^^ "Q . . , x,, . . .) 

(6b) ^^l(. . . , + TT, . . .) = JJ a, ^l(. . . , Xi, . . .) 

where the . . . mean unspecified variables Xi, X2, etc, p = e"^"^ and zi = e~^'^'% i = 1, . . . , L. 

(Note that the factor 11^(^4) to be expected since \E'l only depends on differences 

of spectral parameters. The factor can be absorbed in a redefinition of \E'l, but is 
convenient. The factor 11^(^1) '^j is again expected from the properties of the i?-matrix by 
shift of vr; it could be absorbed in a simultaneous redefinition of the i?-matrix and of \E'l.) 

Similar properties have been observed and (in some cases) proved for models based on 
trigonometric or rational solutions of the Yang-Baxter equation at special points of their 
parameter space [11, 12, 10, 9, 33, 34, 25, 6, 8], except the entries are ordinary polynomials of 
prescribed degree (the main difficulty being to prove this degree). In particular, in the limit 
C — )■ 0, reduces to the eigenvector of the inhomogeneous six- vertex transfer matrix whose 
existence and uniqueness was proved rigorously in [25] and references therein. Therefore, if 
such a solution of (5) exists, it is necessarily unique (for generic p) up to normahzation. In 
principle this normalization might contain a non-trivial function of the Xi, which is why we 
added to the conjecture the fact that the entries have no common factor. So there remains 
only an arbitrary constant in the normalization of "^l, which will be fixed later. 

2.2. Exchange relation. As a direct application of the Yang-Baxter equation, we have the 
following intertwining relation: 

(7) ^l(^| • • • ) Xi, . . .)Ri^i^i{xi^i — Xi) = -Rj^j+i(xj_(_i — Xi)Ti(u\ . . . , Xj, x,_|_i, . . .) 

(see Lemma 1 of [■ ] for the same formula in a similar setting, and its graphical proof). 
Now apply \I'l(xi, . . . , x^) to Eq. (7) and use the eigenvalue equation (5a): 

Tl(m|Xi, . . . , Xj_|_i, Xj, . . . , X2,)-Ri,i4-l(Xj-|_i — Xj)\I'i(. . . , Xi, Xj_|_i, . . .) 

= tL{u\xi, XL)Ri,i+l{Xi+i - Xi)^L(. ..,Xi, Xi+i, . . .) 
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tL{u\xi, . . . ,xl) being invariant by permutation of and F^, commuting with Ri^i+i, 

we conclude by the uniqueness of the solution of (5) that 

Ri^i+i{xi+i - . .,Xi,Xi+i, ■■■) = ri{xi, . . . . .,Xi+i,Xi, . . .) 

where is some scalar function which is a ratio of theta functions, but cannot have a non- 
trivial denominator because it would be a common factor of the "^a, which would contradict 
Conj. 1; so it is a theta function of degree 2 in Xj, Xj+i (and zero in all others, hence a constant) 
with given pseudo-periodicity property; by applying the identity twice and using unitarity 
equation (3), we find rj(xj, Xj+i)rj(a;j+i, Xj) = r{xi — a;j_|_i)r(xj_|_i — Xt). The only theta 
function which divides the right hand side and has the same pseudo-periodicity properties 
as ri{xi, Xj+i) is r(xi+i —Xi); so rj(xj, Xj+i) = ±r(xi+i — Xj). The simplest way to fix the sign 
is to use the C ~^ limit where it is known [25] that the correct sign is +. By continuity in 
C, we have in the end rj(xj,Xj+i) = r(xj+i — Xi), so that 

(8) Ri^i+i{xi+i - Xi)*L(. ..,Xi, Xi+i, ...) = r{xi+i - Xi)*L(. . . , x^+i, Xj, . . .) 

2.3. Spin flip. Next, note that weights a(x) and b{x) (resp. c(x) and d{x)) are exchanged 
by shift of x by vrr. More precisely, we have the following identity: 

Roi{x + vrr) = -p-^ z FiRoiix)Fi 

where z = e~^'^ and Fi is the operator that flips the second spin (of course the same would 
be true with Fi replaced with Fq, since the R matrix commutes with FqFi). 

Applying this to the transfer matrix, we find: 

(9) Tl{. . . , Xi + 7rr, . . = -p~^Zi FiTii- . . , Xj, . . .) 

where we recall that F^ flips spin i. 

As in the previous section, apply \1'l(xi, . . . , x^) to Eq. (9) and use the eigenvalue equation 
(5a): 

Tl{. . . , Xi + vrr, . . .)Fi^! l{- . . , Xi, . . .) = -p'^Zi tL{u\ . . . , Xj, . . .)Fi^! l{. . . , Xj, . . .) 

We have tL{u\ . . . ,Xi + tit, ■ ■ ■) = —p~^ ZitL{u\ . . . , Xj, . . .), as should be, and F^ and Fj com- 
mute, so we conclude as before that 

(10) Fi^L{. . . , Xi, . . .) = /i(xi, . . . , Xl)^l(. . . , Xi + vrr, . . .) 

where /i(xi, . . . , xl) is a scalar function with the following properties: it is a ratio of theta 
functions, but cannot have a non-trivial denominator because it would be a common factor 
of the \I/q,, which would contradict Conj. 1; so it is a holomorphic function, with pseudo- 
periodicity properties determined by shifting one of the Xj by vr, vrr in Eq. (10) and comparing 
with (6); we find 

/i(. . . ,Xi + vr, . . .) = /i(. . . ,Xi, . . .) 

/i(. ..,Xj + vr, ...) = -/i(. ..,Xj,...) j 
/i(. . . , Xi + 2vrr, . . .) = p^"/i(. . . , Xi, . . .) 
/i(. . . ,Xj + 2vrr, . . .) = p~Vi(- --.Xj,...) j 

This fixes it to be /i(xi, . . . , x^) = ce^™^'~'^-'(^') . By rotational invariance, the con- 
stant c is independent of i. Iterating Eq. (10) results in /i(. . . , Xi, . . .)/i(. . . , Xi + vrr, . . .) = 
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p'in^'imxi 2i J2j(^i)^]^ which imposes that c = ±p". In order to fix the sign, we use the in- 
variance by shift of all the spectral parameters and the fact that -F*\E'i = (— with 
F* = Y[i=i Fi to conclude that c'^p~'^(^^i)/2 = i and therefore c = (— p)". 
We finally obtain: 

(11) 



2.4. Wheel condition and recurrence relations. We are now interested in the situation 
where two successive spectral parameters have difference 2?]. In this paragraph, we denote 
to simplify = Tl{u\ . . . ,x,x + 2r], . . .) and = Tl{u\ . . . ,x + 2r], x, . . .) where the two 
specialized spectral parameters are at sites i,i + 1. Applying the intertwining relation (7) 
with Xi = X + 2?7, Xj+i = X, we find: 

A direct calculation shows that R{—2r]) = 2d/^{2ri,p'^)di{2ri^p'^)di{{],p^)P where P is the 

/o o\ 

projector P = 1(1 — P) = 1 j o J^i o ) • Therefore the equality above says that leaves 

Vo 0/ 

ImPj^j+i stable (and that restricted to that subspace it is equal to the projection of T£: 



T; 



L |ImPi,i+i 



' T'" I ^ 
1,1+1-^1, |ImPi,i+i, 



We shall need to check this explicitly. Set s = \ ti) — | it) = -i J be a generator of 
the image of projector P, and compute Ro^i{x)Ro,i+i{x + 2ri)si^i+i, e.g.. 



(-^ I RQ^i{x)Ro,i+i{x + 2r])si^i+i \ ^) 



= {d{x)a{x + 2ri) — b{x)d{x + 21])) \ H) 
+ (a(x)c(x + 2r]) - c{x)b{x + 2r])) \ tt) 
= 



(-^ I Ro,i{x)Ro^i+i{x + 2r])si^i+i \ 



<- 



= (a(x)6(x + 2?]) - c{x)c{x + 2?])) 1 1|) 



8 P. ZINN- JUSTIN 

+ {d{x)d{x + 2r]) - b{x)a{x + 27])) \ l^) 
= r{x + Qri)r{x + 2r])si^i+i 

and similarly with all arrows reversed. Thus, -Ro,i(a;)-Ro,j+i(a; + 2?7)si,i+i = r{x + 6ri)r{x + 
2ri)si^i+i ® lo, and therefore, after shift of x ^ x — u, and use of rj = tc/S to get rid of the 
61], we find 

(12) T+lii^P^ = r{x - u)r{x + 2r] - m)Tl_2 

where it is understood that Ti_2 acts only on sites distinct from i,i + 1. 

Now apply \1'l_2 (with parameters xj except Xi,Xi+i) tensor Sj^j+i and use eigenvector 
equation (5a): 

T/^'l_2(. . .) ® Si^i+i = r{x - u)r{x + 2r] - u)tL-2{u\ . . .)^l-2(- • •) ® s^^^+i 

By definition, tL{u\ . . . ,x, x + 2ri, ■ ■ ■) = r{x — u)r{x + 2r] — u)tL^2iu\ . . .). Also, from Eq. (5b), 
F**L_2 ® = ((-l)"-i^L_2) ® (-Si,i+i) = (-l)"^L-2 ® By uniqueness of the 

solution of (5), we conclude that 

(13) ^Li- ..,x,x + 2r],...) = iJi{x; . . .)^L-2i- ■ ■) ® Su+i 

where, by the same kind of argument as in previous sections, tpi is a theta function of its 
arguments of degree 1 in the Xj, j 7^ z, i + 1 and of degree An in x. 

In order to fix the function ipi, we shall need the so-called wheel condition vanishing 
relation. Let us first consider a special case of it: suppose three successive spectral parameters 
Xi, Xj+i, Xj+2 are of the form x,x + 2r], x + ir]. Then according to Eq. (13), 

Pi,i+i^L(. ..,x,x + 2r],x + 4:r],...)= Pi+i,i+2^L(- . . ,x,x + 2r],x + Ar], . . .) 

= -\1'l(. . . ,x,x + 2r],x + At], . . .) 

But the action of the symmetric group ^3 on ® does not possess the sign 
representation ClS db sub-representation; therefore 

li- ■ ■ , X, X + 2r], X + At], . . .) = 

Now assume all other parameters Xj, j ^ i,i + l,i + 2, are generic; then according to Eq. (3), 
R{xj — Xk) {j i,i + l,i + 2, k = i,i + l,i + 2) is an invertible operator. Applying repeatedly 
the exchange relation (8) to the equality above, we conclude that 

(14) ^l(. ..,x,...,x + 2r],...,x + 4:ri,...) = 

where the location of the three arguments is now arbitrary, as long as the cyclic order is 
respected. This is the general wheel condition, (the equality is true for generic xj, therefore 
for all Xj). 

Finally, using pseudo-periodicity relations (6), as well as flip relation (11), we conclude 
that the wheel condition vanishing relation (14) is valid provided the triplet of spectral 
parameters forms a wheel x,x + 2r], x + At] modulo tt, ttt (not just 2ttt\ a crucial technical 
point which will be used repeatedly below). 

We can now come back to our recurrence relation (13). On the left hand side, we notice that 
as soon as one of the xj, j + is equal to x — 2?7 (mod vr, vrr), a wheel is formed and "^/l 
vanishes. Therefore ipii^'i ■ ■ •) contains factors H^x^i i+i) ''^i(^~2r7 — Xi;p^)'(?4(x — 2?] — xi;p^); 
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moreover these exhaust its degree, and noting that these factors can also be written up to a 
multiphcative constant as 'di{x — 2rj — Xi;p), we can rewrite Eq. (13) 

(15) ^L(...,x,a; + 2r/,...) = est JJ ??i(a; - 2r/ - x^; p) ^l_2(. • •) ® ^^M+i 
More exphcitly, it means that 



^L;ai,...,«L \xi+i=Xi+2ri 



ai = tti+i 

CStaiY[j(j^i^i+l)1^l{Xi -2r] - Xj;p) ^L-2;ai,...,ai_i,ai+2,...,aL 7^ 

The constant remains undetermined at this stage, since we have not fixed the normahzation 
of yet. 

A similar recurrence relation can be written for Xj+i = Xi + 21] + ttt (the non-zero result 
occurring when a, = aj+i), but we shall not need it. 



3. Partition function 



In the rest of this paper, we denote ^{x) := {}i{x;p) and i!^k{x) := 'i9fc(x;p), k = 2,3,4. 
Since the contents of this section are not expected to generalize outside r] = n/S, we shall 
use 31] = (mod vr) to replace 2// with —rj whenever possible. 

3.1. Definition. We now introduce a quantity that naturally generalizes the squared norm 
of the XYZ ground state (Eq. (4)) to the inhomogeneous case: 

Zi(xi, ...,xl) = (^L(-a;i, • • • , -xl) \ ^L(a;i, . • . , xl)) 

where we have used the (real) scalar product: ($ 1$') = X]ag{t4,}^ ^a^'a- 

Zl{xi, . . . ,xl) has the following interpretation: it is the "partition function" of the inho- 
mogeneous eight-vertex model on an infinite cylinder. Indeed, assuming that we are in a 
regime of parameters where is associated to the largest eigenvalue of the transfer matrix, 
\E'i(a;i, . . . , Xl) corresponds to the partition function on a half-infinite cylinder (pointing up- 
wards) with given arrows at the boundary at the bottom. A vertical mirror symmetry of 
the eight vertices correspond in the weights (1) to x — )■ —2i] — x and a change of sign of the 
weights a and b, the latter being irrelevant with periodic boundary conditions. So the parti- 
tion function of the other half-infinite cylinder (pointing downwards) is \E'l(— Xi, . . . , xl) {"^l 
only depends on the differences of its arguments so the —2i] term is irrelevant). We mean 
partition function in the following sense: a one-point correlation function will be expressed 

as {O) = \ Ea,/3e{t,4}^ ^L;a(-2;i, . . . , -Xl)0 ^^fi^ L-fiix^, ...,xl). 

3.2. Pseudo-periodicity. According to its definition and (6b), is invariant by Xj — )■ 
Xj + vr for any given i. Furthermore, 

Zl(. . . ,Xi + vrr, . . .) = ^ (^j:(. . . , -Xj - vrr, . . .) | ^'i(. . . , Xj vrr, . . .)) 

ae{t,;}^ 

= (-p)"e'™(-"^-"")+'^:'(^')"^" (^i(. . . , -Xi, . . .) I 

(_p)-ne-2™-.+iE,(^.)X.jr. I vl/^(. . . ,0;,, . . .)) by Eq. (11) 
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where it is reminded that zj = e . 

We reach the conclusion that Z^, is a theta function of degree 2n and nome p (as opposed 
to for "^l) in each variable Xi. 

3.3. Symmetry. Given i = 1, . . . , L — 1, we can use the exchange relation (8) and unitarity 
relation (3) to write 

Zl{xi, . . . 

= . . . , -Xi, . . . , -Xl) I ^LiXi, Xi+i,Xi, . . . , Xl)) 



riXi — 1 



, . . . , Xi, 



Ri,i+l{Xi — Xj+i) 



-^l{xi, . . .,Xi,Xi+i, ...,Xl) 



• • • 5 ■"^JlJ 



'"(^^i+l ~ Xi) 



'^L{Xi,...,Xi,Xi+i,...,XL)) 



Zl{xi, ...,Xl) 



where in the intermediate step we also used the fact that the R matrix is self- adjoint. We 
conclude from this calculation that is a symmetric function of its arguments. 

3.4. Recurrence relation. The recurrence relation (15) for implies one for Z^: 
..,x,x + 2r]) = (^i(. . . , X, X + 2r/) I ^i(. . . , -x, -X - 2r/)) 



L-2 



cx Y[^x-2r]-Xi) (^L 

i=l 
L~2 

CX JJ 79(x - 2r/ - X,) (^L 



SL-1,L I *l(- ..,-X,-X- 27])) 

sl-i,l I Pl-i,l I ^l(- . . , -X, -X - 2?7)) 



L-2 



cx JJ 7?(X - 2r/ - Xi) (^'l_2(. . .) ® Sl-1,L | *l(. ..,-x-27], -x)) 

by Eq. (8) with Xl_i = — x, x^, = — x — 2?] 

cx Y\ — 27] — X.i)§{{ — X — 27]) ~ 27] — { — Xi)) 

{^L-2{Xl, . . . , XL-2) ® Sl-1,L I '^L-2{-Xl, -XL-2) ® Sl-1,l) 



i=l 



L-2 



i=l 
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L-2 

OC 

i=l 

where oc means equal up to a multiplicative constant. At this stage, we fix the normalization 
of in such a way that this constant disappears in the recurrence formula for Zl, which 
becomes: 

L-2 

(16) Zl{. . . , X, X + rj) = Yl'i^'^ix - V - Xi) Zl-2{- ■ •) 

1=1 

where we have also shifted x — > x + 77 and used 3?7 = (mod vr). 

Combined with the symmetry in its arguments, the recurrence relation (16) satisfied by Zl 
means that we can express its specialization at Xi = X2 ^rj, ... ,xl ±7] in terms of Zl_2. So 
we possess 4n values of as a function of Xi, since it is a theta function of degree 2n, these 
relations are more than enough to determine Z^ inductively (say, by Lagrange interpolation). 

3.5. Half-specialization. At the moment, we do not know how to solve in a closed form the 
recurrence relation above. However, note that we have twice as many recurrence relations 
as needed to determine Zl- This suggests to "half-specialize" Zl in such a way that the 
number of recurrence relations now matches the degree. 

Explicitly, assume Xi+„ = — Xj, i = 1, . . . ,n, and xl = 0. After such a specialization, Zl 
is an even function of xi, . . . , x„, and it has a double zero at Xj = ±77, i = 1, . . . ,n. Let us 
check the latter statement carefully. Since Zl is a symmetric function of its arguments, let 
us assume that we order them as x, 0, — x . . . (where x is one of the Xj) and that we send x 
to — 2?7 (which is equal to 1] modulo vr). Then it is clear that \E'l(x, 0, — x, . . .) forms a wheel 
and therefore Zl{x, 0, — x, . . .) vanishes. However, \E'i(— x, 0, x, . . .) does not vanish, so that 
to show that the zero is double, we need to go further. Write 



d 

(17) — Zl(x,0, -X, . . .)|^=_2r, 



— ^'l(x, 0, 2r/, . . .) |,.=_2^ - — ^L(-2r/, 0, x, . . .) U=2„ 



^i(2r/,0,-2r/,...) 



Now apply recurrence relation (15) to \l'/,(x, 0, 27], . . .); we find that it is proportional to 
some vector at sites j 7^ 2,3 tensor ^2,3, and therefore the same its true of its derivative 
w.r.t. X (one can be more explicit using x = —2r] but we shall not need it). Similarly, 
^'i;,(— 27], 0, X, . . .) and its derivative w.r.t. x are equal to Si^2 tensor some vector at other 
sites. 

On the other hand, applying the exchange relation (8) to \E'j;,(2?7, 0, —2?]) at i = 1 implies 
that Pi^2^l(277, 0, —2rj, . . .) oc \I'l(0, 2?7, —2rj, . . .) = since a wheel is formed. Similarly, the 
exchange relation aX i = 2 implies that P2,3^l(2?7, 0, —2rj, . . .) = 0. 

We conclude that the expression (17) is zero by inserting P2,3 (resp. ^1,2) in the first (resp. 
second) term. Therefore, taking into account evenness, we can write 

n 

(18) Zl(xi, . . . ,X„, -Xl,. . . , -X„,0) = JJ^9^(Xi - 7])^'^{Xi+7])y.n{Xl, . . . , X„) 
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where Xn{xi, . . . , Xn) has the following properties, as a direct consequence of the correspond- 
ing properties for Z^,: 

• X„ is a symmetric function of its arguments, and an even theta function of degree 
2(2n — 1) in each. 

• It satisfies the recurrence relations: 

n-2 

(19) X„(. ..,X,X + 7])= (p\x)(p\x + 7])Yl -Tj- Xi)l^\x -7] + Xi) X„_2(. . .) 



1=1 

n-1 



(20) X„(. . . , /3fe) = n ^--i{- ■■) A; = 2, 3, 4 



i=l 



where ip{x) = = ^'&2{x)'&3{x)'&4{x) , n = ^^^^^^ = i92(?7)i?3(?7)i?4(?7) and H2 = 

1, ^3 = ^4 = — i/^""^ (i/ = e~^'^*/^/y^) are pseudo-periodicity constants. 

where /32,3,4 are representatives of the three solutions of 2/3^ -|- 77 = (mod tt, ttt) excluding 
77, namely, ^2 = 7r/2 77, ^3 = 7r/2 'kt/2 77, ^4 = 7rT/2 77. 

We now have at our disposal the specializations xi — ±Xi±r], ±/3k, i = 2, . . . , 7i, A; = 2, 3, 4, 
that is 4(71 — 1) -1-6 = 2(271-1-1). An even theta function of degree 2(27i — 1) being determined 
by 2 X 271 values, we have enough recurrence relations to determine X„. 

3.6. Solution as PfafRans. We first introduce the function: 

which has the following properties: 

• It is symmetric function of x, y, and is an even theta function of degree 2 in each. 

• It satisfies the following recurrence relations: 

(21) A2(x, x + 7]) = -v^ d2.{x)^-i{x + r])^4.{x)di{x + 77) 

(22) k2{x,l3u) = pHl{ri)dl{x) A; = 3, 4 
Next we claim the following: define 

(23, "-'--'-")-,n„ .(Jt];^^+., -^^- 

where 

h(a;, y) = 79(77 + x — y)^{r] + x + y)-^{r] — x — y)^{r] — x + y) 
and M„ is a skew-symmetric 2771 x 277i matrix, m — [ri/2], given by 

i{xi, Xj) n even or i,j < 2m 

-1 n odd, i — 2m, j < 2m 

1 n odd, j = 2m, i < 2m 

71 odd, i — j — 2m 



(24) (M,) 
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and 

(25) f(x, y) = ^i^-yn^ + y)Mx^y) 

h[x, y) 

Also define 

(26) Bn{Xi, ...,Xn)= A„+i(xi, ...,Xn, /32) 

Then: 

• An (resp. B„) is a symmetric function of its arguments, and an even theta function 
of degree 2(n — 1) (resp. 2n) in each. 

• They satisfy the recurrence relations: 

(27) A„(. . .,X,X + T]) = -Z/^ ■&3{x)l93{x + ri)-&4,{x)-&4,{x + 7]) 

n-2 

Y\_'&'^ix - 1] - Xi)'d^{x -Tj + Xi) A„_2(. . .) 
1=1 

(28) B„(. . . , X, X + r]) = -z/^ ^l{x)^l{x + 27])^-i{x)^:!.{x + 7])^i{x)^A{x + t]) 

n-2 



^^{P{X -T] - Xi)'&'^{x - Tj + Xi) B„-2(- • •) 



1=1 

(29) A„(...,/32)=B„_i(...) 

n-l 

Mf^ \ A 

i-n-1 



(30) B„(. . . , /32) = -^^^l{r]¥l{v) n ^2(^.) A. 

i=l 

n—1 

(31) A^. . . , /3,) = v^dM{^H'n)Y-'^^ n ^^(^^) • •) ^ = 3' 4 



i=l 



(32) B„(. ..,/?,) = z/"+it92(r])(z/^?2(r]))-(-^)" !l!pM TT ^2 _ _) ^^3^4 

where conventionally Aq = Bq = 1. 

Let us show for example (27) for n = 2m even. Assume that X2m~i = x approaches X2m + 
7] = x' . Then the matrix element (M„)2m-i,2m develops a pole: (M„)2m-i,2m oc l/(x — x'); 
and the other entries (M„)jj (i < j) remaining finite, the only relevant contributions to the 
Pfafiian are those pairing 2m — 1 and 2m, so we immediately have 

A / I N TT h(xj, x)h(xj, X + r]) , na / \ 

An • . . , X, X + r/ = M — — -— ■ — — ■ ■ — -A2 X, X + r/ A„_2 . . . 

v\Xi — x)i}[Xi — X — ri)V[Xi + x)iy[Xi + x + rj) 

where we have cancelled all factors in common to A„ and A„_2- 

Now the remarkable phenomenon (using in a crucial way t] = vr/3) is that there are 
compensations in the product, which simplifies to YYi=i '&'^{x — i] — Xi)'d'^ {x — 1] + Xi). Finally, 
we use Eq. (21) for A2 to reproduce the remaining prefactors on the r.h.s. of Eq. (27). 

The other equations follow from similar reasonings. 
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Finally, we find that A„B„ satisfies all the recurrence relations of X„, or more precisely, 

X„(X1, ...,Xn) = (-Z/^/«^)""A„(xi, . . . , Xn)Bn{Xi, . . . , X„) 

3.7. Further factorization as determinants. Consider the following elliptic version of 
Tsuchiya's determinant [30, 20]: (see a similar determinant in [14]) 

m 2m 



n n H^i'^j) 



(33) H2m(2^1; • • • 5 ^m] ^m+l^ • • • j X2m) — ' dct 



n ^{x. - x,)d{x, + X,) J^f^j^,^ h(x„ X,) 

l<j<j<m 
or 

m+l<j<j<2m 

Conventionally, Hq = 1. Note that = 1 as well. 

The expression of has the disadvantage that it is only (apparently) symmetric in the 
variables {xi, . . . , x^} and {x^+i, • • • , X2m} (we shall show below that, thanks to rj = 7r/3, it 
is in fact symmetric in all variables). In terms of each, it is an even theta function of degree 
2(m - 1). 

It is not hard to see that satisfies the following recurrence relation: 

(34) Y{2m{- ■ ■ ,x;. . . ,x + 7]) = JJ ?9(x - ?7 - Xi)^?(x - ?7 + Xi)H2m_2 (...;.. .) 

i=l,...,m— 1, 
m+l,...,2m 

and the same if one exchanges x and x + rj. These are the usual recurrence relations satisfied 
by such determinants, as in the classical case of the Izergin-Korepin determinant [18, 16, 17], 
and similarly to the Pfaffians of Sect. 3.6; the reasoning to derive Eq. (34) is identical - a 
pole develops in one of the entries of the determinant, reducing it to a determinant one size 
smaller. 

Now consider the function 

^2miXl^ • • • ) X2m} ~ H2m(3^1; • • • ; ^m; ^^m+l; • • • 5 2;2m)H2m+2 (^^1 ? • • • ? Xjyi, (3^^ Xm,+1, • • • ; X2rm P4} 

is an even theta function of its arguments, of degree 2(2m — 1), and using Eq. (34), it 
satisfies the same recurrence relation (27) as A2m, for say Xi = ±Xj ± r/, j = m + 1, . . . , 2m. 
Thus the function is known at 2 x 2m values of Xi, which determines it uniquely. Combined 
with A2 = v42 = 1, we conclude by induction that Ag^ = A2m- 

An interesting corollary is that B.2m is a symmetric function of all Xi, . . . , X2rrt together. 

Similar arguments can be made for A2m-i and B„. Together, we find 

(35) A2m(a;i, . . . , X2m) = ii2m{Xi, . . . , X2m)H2m+2 (a^l , . . . , X2m5 /^Sj (^i) 

(36) A2^ _l(Xi, . . . , X2m-l) — ii2m{Xi, . . . , X2m-1, f33)^2m{ 

Xi , . . . , X2rrt— 1 ) 

(37) B2m{xi,... ; X2m ) — H2m+2( 

Xi , . . . , X2m; 1^2, /33)H2m+2( 

Xi , . . . , X2m; 

(38) B2^_i(xi,... ; 3;2m--l) — B.2m{Xi, . . . , X2m-1) /32)H2m+2 (^^l , . . . , X2m-1, /32; /^Sj (^^4) 

These are the only 8 possible specializations at /32,3,4, corresponding to subsets of {f32, /^s, /34}, 
since applying any such specialization twice amounts to the shift n — > n — 2. 
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3.8. Alternative determinant formula. Here we follow the same general method as in 
[29] (see also appendix B of [11]). Define 

Hx,y) 

g{x, y) is an odd elliptic function of x and y. One further observes that for all x, y, 

g{x, y) + g{x + ri,y) + g{x + 27], y) = 

and similarly for y. Therefore, detj , g( function of any of its arguments, satisfies 

the same three-term relation. Now define 

S2miXi,...,X2m) = JJ l^{Xi - Xj) JJ + Xj) HamlXi, . . . , Xsm) 

l<i<j<2m l<i<j<2m 
2 

nTT - Xj + kr])-d{xi + Xj + kr]) det g(xi, x,) 
l<i<m 

l<i<i<m k=0 m+l<j<2m 
or 

m+l<i<j<2m 

Since the prefactor is invariant by Xi + r] for any i, we have the same three-term relation 
for S2m- In summary: 

• is a skew-symmetric function of its arguments Xj, and an odd theta function of 
degree 6m in each. 

• It satisfies 

S2m(- . . , X, . . .) + S2m(. . . , a; + ?7, . . .) + S2m(. . . , a; + 2?], . . .) = 

The space of odd theta functions of degree 6m is of dimension 3m, a possible basis being 

s^{x) = e^'^^'d^ik-KT + 6mx,/") - e-'^'^'^^^iik-nT - Qmx,p^"') k = 0,...,3m-l 

Sk{x) satisfies the relation Sk{x) + Sk{x + ri) + Sk{x + 2ri) = iff A; 7^ (mod 3). The sequence 
(1, 2, 4, 5, ... , 3m — 2, 3m — 1) = {ki, . . . , k2m) is of cardinality 2m, which is the number of 
variables of S2m, so we conclude that 82™ is proportional to the "Slater determinant" 

(39) S2m(Xi, . . . ,X2m) OC dct SkAXi) 

i,j=l,...,2m 

We shall not need the proportionality constant. 

3.9. Uniformization. Although the formulae above are simple to derive, they are a bit too 
cumbersome to be used, especially in the homogeneous limit. Since all functions we consider 
are theta functions of definite parity, there is a rational uniformization, and we use from now 
on the following parameterization: 

2\-l/3 l^'^jx) 



"^^""^ ^{X - l])d{x + T]) 

In terms of the original Boltzmann weights (1), we have (1 — C'^)w{x) = ^"'^^(^^^l^^Zli^^ 



Note the special values 



w{/32) = -^ = J2 wiPs) = = J3 w{/34) = = Ja 
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which explains the labelhng we have chosen for the couphng constants </2,3,4- 

This parameterization has the advantage that the wheel condition becomes simple to ex- 
press: three spectral parameters form a "wheel" ±x, ±{x + ri), ±(a; + 2?7) iff the corresponding 
variables w, w', w" satisfy 

3 + 



(40) 



w + w' + w" 



w w'w" 



1-C 
1 



1 

and therefore, two parameters form a "2-string" ±x, ±{x + ri) iff the corresponding variables 
w,w' satisfy h{w,w') = 0, where 

(41) h{w, w') = l-{3 + C^) W + (1 - C^)ww'{w + w') 

This formula allows to rewrite the recurrence formulae in this new parameterization, but 
due to the fact that it is quadratic in w and w', the result is somewhat cumbersome and we 
shall not write it explicitly. 

We also redefine the functions by dividing them by a "reference" even theta function of 
degree 2 to the appropriate power, here ■(9(x — ri)d[x + r/), and absorbing some constants in 
the normalization. That is, we define 



An{w{Xi), . . .,w{Xn)) = a„ 
Bn{w{Xi), . . .,w{Xn)) = 6„- 



Bn (2^1 , . . . , Xn) 



where a„ and bn are constants which are implicitly defined by the expressions below, and 
whose explicit expression we shall not need. 

In particular, 

A2{w, w') = ww' - {w + w') + — 

1 - c 

and if we define 

{w - w')A2{w,w') 

f{w, w = 

h[w, w') 

which is such that f{w{x),w{y)) = —(q/y/p i{x,y), then we have: 
(42) A4w,,...,wn)= n ^^^^^ PfM„ 

l<i<j<n * J 

where M„ is identical to M„, except entries f{xi, Xj) are replaced with entries f(wi, Wj); and 

Bn{Wi, . . . , Wn) = An+l{Wi, ...,Wn, J2) 

as well as 

Xn{Wi, ...,Wn) = 2''+M„(wi, . . . , Wn)Bn{Wi, . . . , W„) 

where the numerical coefficient has been adjusted so that in the rational limit, the normal- 
ization of Xn coincides with the one discussed in Sect. 1.3. 
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A further advantage of this new normahzation is that An and -B„ are polynomials in 
Wi, . . . ,Wn and also of (, up to a conventional denominator in powers of 1 — which we 
have added for convenience. 

Similarly, we can define 

m 2m 



n n 



, i=l j=m+l 
H2m{Wl, W2m> = dct 



l<i<j<m 
or 

m+l<i<j <2m 



and then the relations (35-38) expressing A, B in terms of H remain the same; more com- 
pactly, one can write: 

SC{J2,J3,J4} 
\S\=n (mod 2) 

There are various alternative formulae, for example 

m 2m 

n n ^(^i'^i) 

TT r 7 TN i=lj=m+l A2{Wi,Wj) 
H2m+2{Wl, . . . , W2m, J3, Ja) = dct — 

n (^^-^^)mim<2m ^(^"^^^ 

l<i<j<m 
or 

m+l<i<j <2m 



Finally, the transformations 



C-1 

generate the group of permutations of the three coupling constants J2, J3, J^. Via the uni- 
formization w{x) = W!^{x), this translates into the symmetry of permutations of non-trivial 
solutions of 2x + 1] = 0. The function h{w, w') = h(^{w, w') itself possesses this symmetry, in 
the sense that 

h^(w, w') = h^c{w, w') h^{w, w') = /i(^+3)/(^_i) {-^—w, -^—w') 
which is consistent with w^{(33^i) = W-,j(/34,3) and 

W(c+3)/ic-i){f32) = -1/2 = •^4 

(R^ C-1 C-1, 

W(C+3)/(C~l)lP3) = 2{( + l) ^ ^ 
W^(C+3)/(C~l)lP4j = = ^—-^2 
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4. Homogeneous limit of the partition function 

4.1. Summary. The homogeneous hmit is obtained by setting all spectral parameters equal; 
in the half- specialized partition function X„, this is achieved by sending all to zero. 

In this section, we use the following notation: we omit parameters that are set to zero, 
e.g., if 2m = -^2m( p, ■ -j ,0) - This is unambiguous because the total number of variables is 

2m 

given in subscript. Here are some values of for ^ = 0, 1, 2, 3: 



H2m = 1; 


1, 


3+c^ 


26 + 29C^ + 8C^ + 






2™ ^H2miJ2) = 


1, 




143 + 99C^ + 13C^ + C*^ 






H2m{Jz) = 


1, 


2 + C + C', 


11 + 12C + 21C^ + 10C^ 


+ 7C' + 


■2C' + C' 


H2m{JA) = 


1, 


2-C + C', 


11- 12C + 21C^- lOC^ 


+ 7C^- 


•2C' + C' 


2™" ^H2m{J2, Js) = 


1, 


5 + 2C + c', 


66 + 63C + 81C^ + 30C^ 


+ 12C'- 


+ 3C^ + C' 


2™" ^H2m{,J2, Ja) = 


1, 


5-2C + C', 


66 - 63C + 81C^ - 30C^ 


+ 12C' 


-3C' + C' 


H2m{J^i Ji) = 


1, 


l+c^ 


3 + 9C' + 3C' + C' 






2"^ ^H2m{J2^ Js, Ja) = 




3+c^ 


21 + 39C^ + 3C^ + 







We recognize the reciprocal polynomials of those occurring in conjecture E of [ ]: the 
correspondence of notations is that for m > 1, H2m = C^'"*'™~^''Q'm--i(l/C)) H2m{J'i) = 

All other sequences can be obtained by permutations of the {J2, J3, J4}, and can therefore 
be obtained by iterating the transformations C — ^ — C ^-^id ( — )■ as explained at the end 
of last section. All the properties listed in conjecture E of [5] can thus be checked on the 

H2m- 

If we recombine the H2m in pairs to form An and we recognize the reciprocal polyno- 
mials of the Sn of [ ]: (see their appendix A) 

(43) A„ = C'LnV4j,^(l/^2) 

(44) Bn = (2/3)X'L("+'^'/'Js-n-i(l/C') n > 
from which we conclude 

X„ = 2^+'AnBn = 2(4/3)"C^""''^s„(l/C')s_„_i(l/C') 

which coincides with the expression given in Conjecture 1 of [5] up to the factor of two (which 
is due to our slightly different way of lifting the two-fold degeneracy: F^'^l = (— 1)"'^l 
effectively duplicates every entry of compared to [ ]). 

See also appendix A for an explanation of the constant terms of the various polynomials 
above. 

In the rest of this section, it is convenient to denote a = 1 — C^. We shall show that the 
various polynomials above satisfy (differential) recurrence relations. 



SUM RULE FOR THE EIGHT- VERTEX MODEL AT r? = 7r/3 19 

4.2. Linear relations. We first derive certain linear relations satisfied by H2m{wi, . . . , W2m)- 
We shall need them to relate the various derivatives of at = 0. 

Define 

(45) D2m{yJi,...,Wm;Wm+i,---,W2m)= det g{Wi,Wj) 

l<j<m 
m+l<j<2m 

with g{u, v) = j^^', we recall that h{u, v) = 1 + uv{a{u + f + 1) — 4). In other words, 

m 2m 



n n ^(^i'^i' 



(46) H2m{Wl, . . . , W2m) = -D2m(ti'l, • • • , Wm] Wm+1, • • • , W2m) 



l<i<j<m 
or 

m+l<«<i<2m 



Also define 



(47) 5*2^(^1, • • • , U;2m) = Yl {U^i- Wj)H2m{Wl, . . . , W2m) 

l<i<j<2m 

4.2.1. A first order differential/ divided difference equation. We start from the following iden- 
tity, which can be checked directly: 

p{u)dug{u, v) + p{v)d,,g{u, v) 

+ 2(1 - a){8 + a)dag{u,v) + {(t{u) + a{v))g{u,v) + (5„ + 5/)g{u,v) = 

where p{u) = (1 + 2m) (4 — 6u + ua + u^a), a{u) = 5u{a — 4 + Aau), du is the usual partial 
derivative and 5u is the divided difference operator: 6u4>{u) = for any function 

Then, one can easily prove starting from (45) (for example by writing D2m as a sum over 
permutations and grouping together the summands for values of the index connected by the 
permutation) 

2m 



(48) Yl + ^(^i) + O + 2(1 - «)(8 + «)9a 



In principle, by using relation (46), one can reformulate this identity in terms of H2m, but 
the result is not particularly illuminating and we shall not need it. 



4.2.2. A second order differential equation. Starting from the differential equation satisfied 



by 1^3, namely, (|4 + 4p|-)793(x; p) = 0, we find 



92 _ d 



dx"^ dp 
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According to Eq. (39), this implies that 

(49) ( ^ ^ + 24mp— j S2m{Xi,..., X2m) = m{Qm^ - l)S2.m(2;i, • • • , X2m) 



After switching to our rational parameterization and from p to a, we find the following 
equation for S2m- 

(2m \ 
^{j2{wi)di + -fi{wi)di + -fo{wi)) + 24ma{l - a){8 + a)da j S2m{wi, . . . ,W2m) =0 

where the coefficients have rather unpleasant expressions: 
7o(m;) = w'^{3Qa^ + h^o?m^ - 90a2m) + w{ma^ - 144a + 12a^ir? - 2%%air? - \^2Q?m 

+408am) + 64m^ - Vla^m + 96am - 192m + 80 + ho? - 40a + lOa^m^ - 20am^ 
7i(w) = 6(w'^(6a^ — Ga'^m) + w^(9a^ — 36a — lOa^m + 40am) 

+w(3a^ - 24a - 4a^m + 12am - 32m + 48) - 6a) 
72 (u^) = 6w{aw — 4)(a + aw"^ + 2aw — Aw) 

4.2.3. Homogeneous limit. We now take the homogeneous limit in two steps: we first send 
wi, . . . ,Wm to u and w^+i, • • • , W2m to v and then expand around u,v — 0. H2mi being a 
symmetric function of ifi, . . . ,W2mi only has one independent first derivative (resp. two in- 

9 n 



dependent second derivatives), which with our specialization correspond to ^-f^2m = ^-^^2 

(resp. -§[pH2m = -§;pH2m and g^i^2m)- 

Taking this limit in Eqs. (48) and (50) is a rather tedious procedure which we shall not 
describe in detail. Expanding to first non-trivial order these equations produces the same 
result, namely the first equation below. This equation is a first order differential equation, 
and so we can differentiate it once w.r.t. a, resulting in a second order equation (second 
equation below). Expanding to the next order Eqs. (48) and (50) produces two distinct 
second order differential equations. Finally, we find: 



(51) 






2(2m+l) 



2{4m+l) 

2m(l-a)(a + 8) 2(4m + l) 



\2 (4m +m — 2)a — 2m(4m + l)a 



2m(l-a)(a + 8) m(m^ -m+1) (2 + a) 



m(2m + l)(a-4)'^ 



2m(l-c«)(c« + 8) 
-an 





(m-l)(2 + a) 
(m— 1) 



m^(m — 1)(4 — a) 



-(m-l)m'^ (2m + l)(c« — 4)c« 



rH2„ 
tH2„ 
:H2. 



dcx.- 



\ 



S 



da 



H2r, 



H2„ 



There are 4 relations for seven derivatives, so they can all be expressed in terms of deriva- 
tives w.r.t. a only. 
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A similar reasoning can be made when all variables are specialized to except one, or two, 
or three, are specialized to a subset of {J2, J3, Ji}- The result is given in appendix C. 

4.3. Bilinar recurrence relations. We now show how to derive differential bilinear recur- 
rence relations for and its variants. In fact these relations were mentioned, but not 
written explicitly, in paragraph 3 of [■]. 

Similarly to the previous paragraph, we first consider the quantity 

_ 2 / 1 d'+^ 

H2m{u,...,u,v,...,v) = g{u,v) det ( ., ., ^ . 9{u,v] 

m m 

A standard application of the Jacobi-Desnanot identity (see [!!)] for the simpler case of 
the Izergin-Korepin determinant) to the determinant in the right hand side produces the 
Toda lattice equation: 

m+i) 

(m, . . . ,f, . . .)H2 (m— 1) [U^ . . . ,V, . . 

■\ogH2m{u,. . .,v,. . .) = -——\ogg{u,v) + 



rn? dudv ' ' ' dudv ' g{u, vYH2m{u, . . . , f , . . .)2 

The left hand side involves first and second derivative of i^2m(w, . . . ,f, . . .), which at m = 
V = Q can be reexpressed in terms of derivatives w.r.t. a thanks to Eq. (51). The result is: 

(52) C^H2 (m+l) H2 (m—1 

where all derivatives are w.r.t. a, and 
Co = 4a(4m - l)(4m + l)2(4m + 3) 
Ci = 4(a - ifaia + 8f{Am + 1)^ 
C2 = 4(a - ifaia + 8f{Am - l)(4m + 3) 

C3 = 2{a - l)(a + 8){a^ + 28a + 2A.c?ir? + 304am^ - 256m^ + 20a^m + 208am - 192m - 32) 
C4 = ^o? - 24a + 4a^m^ - lOOSaW + 3408am^ + 512m^ - 4a^m^ - 984a2m^ + 4032am^ 
-128m^ + o?n? - U2a^m^ + 1028am^ - 320m^ - a^m + 28a^m - 44am - 64m 



Note that contrary to Eq. (51), Eq. (52) is a closed relation allowing to compute inductively 
the H2m as polynomials of a = 1 — C^. 

A similar computation produces differential recurrence relations of the same form for the 
other factors of X„. The coefficients are given in appendix C. Together, they allow to 
compute the full squared norm X„ inductively. 

5. Conclusion and prospects 

In this paper, we have considered the inhomogeneous eight-vertex model with periodic 
boundary conditions in odd size and crossing parameter 77 = 7r/3. We have provided a ba- 
sic setup for the computation of the fully inhomogeneous generalization of the ground state 
eigenvector of the XYZ spin chain, and then went on to compute the partition function on an 
infinite cylinder, which generalizes the squared norm of the ground state eigenvector, when 
the spectral parameters are "half-specialized", i.e., form pairs x,—x. We have provided a 
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variety of explicit expressions for this partition function in terms of Pfaffians and determi- 
nans. Interestingly, one can then obtain self-contained expressions in the homogeneous limit 
for the squared norm, without any more reference to the inhomogeneous case, by allowing 
differentiation w.r.t. the variable parameterizing the line t] = ir/S (elliptic nome, or (). These 
expressions take the form of bilinear differential recurrence relations (cf Eq. (52)). 

In order to derive such differential relations, we have used certain differential (and divided 
difference) relations satisfied by the inhomogeneous partition function. In fact, we have 
strictly limited ourselves to the relations that were needed for our purposes, but it seems 
that this is only the tip of the iceberg: one should investigate in more detail the structure of 
the set of such equations. It would be interesting to understand the role of the full symmetry 
of arguments of the Izergin-Korepin type determinant (33). 

Note that we have not been able to obtain an expression for the fully inhomogeneous 
partition function, but if we compare to the work of Rosengren for the 8VS0S model [ ] 
there is also no simple expression for the fully inhomogeneous partition function. Inversely, 
it would be interesting to see if the "half-specialization" trick helps in this context. More 
generally, as noted in [")], there are many ressemblances between the work [■2C] and our 
present setup, which should be clarified. 

It is clear that the present methods should allow to compute more quantities such as 
individual entries of the ground state, or certain correlation functions (see the recent work 
[7] in the XXZ setting). 

Some more directions which should be explored are: the relation to the quantum Knizhnik- 
Zamolodchikov-Bernard (gKZB) equation and to the gKZB heat equation [ ], which should 
be the right framework for part of section 2, especially in view of a generalization to arbitrary 
rj; the connection to nonsymmetric elliptic Macdonald polynomials; the use of matrix model 
techniques to analyze the determinants of Izergin-Korepin type found here, as in [32]; and 
the meaning of the connection to the Painleve VI equation, which is emphasized in [4, 5]. 

Finally, it would be interesting to find a combinatorial interpretation for the (positive 
integer) entries of the polynomials of section 4.1, beyond their value at ( = 0. 

Appendix A. The C trigonometric limit 

The trigonometric limit is obtained by sending ( to 0. The Boltzmann weights (1) of the 
eight- vertex model turn into those of the six- vertex (the weight d go to zero). In this limit 
the results of this paper should be closely related to the computations of [!•")]. Note that 
the "quadratic" sum rule considered here was actually not computed in [15] - instead the 
quantity \1/ 0,(2:1, ... , zl)"^ was used there. However, the same argument of degeneracy of 
the scalar product allows to conclude that 

(53) Zi(xi, ...,xl) = 3-"(^ ^aixi, . . . , XL))iJ2 ^"(-^1' • • • ' -^l)) 

a a 

= 3"''%y^(2;i, . . . ,2;i)sy^(zf\ . . .,zl^) 

where sx is the Schur function associated to partition A, and = ([(L — 'j)/2j)j=i j.- 
In the homogeneous limit, sy„(l,...,l) = 3»("-i)/2 YYj^i {2j)\{2i-iy. ^^'^ together we have 
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= Aht{L), where Aht{L) = 1,3,25,588... is the number of Half- Turn Symmetric 
Alternating Sign Matrices [20, 23]. 

The half-specialization of section 3.5 produces the following factorization: 

n 

(54) Sy^(l, Zi, Z^^, ...,Zn, Z''^) = + Zi + Z^^) XyAzU • • • , ^n)Xy„+i (^^l , . . . , Zn, W) 



i=l 

where xa is the symplectic character, defined by: x\{zi, ■ ■ ■ , Zn) = — ' n-j+i — -n+j-is — , 

aet{z^ —z^ ) 

and oj = e*'^/^; this formula can be proved by induction, or can be seen as a byproduct of 
this paper, as we now show. 

In the limit — )■ 0, the parameterization w is related to the multiplicative spectral param- 
eter z by w = {z — 1)^/ (1 + 2; + z^); this way we find 

_ 9{z^ + zz' + z'^){l + zz' + z^z'^) 

{i + z + z^Yil + z' + z'^Y 

The denominator factors out of Pfaffians and determinants. 

A.l. PfafRans. We now recognize the Pfaffian An (Eq. (42)) in even size: 

At ^ Q^TT TT ^{z^i+ZiZ^ + z^^){l + ZiZj+zlz^^) 



i=l l<i<j<2m 



{1 + Zi + Zf){l + Zj + zJ){Zi - Zj){l - ZiZj) 

pf {Zi ~ Zj){l - ZiZj) 



{zf + ZiZj + z]){l + ZiZj + zfz^) 

which up to some prefactors is exactly the Pfaffian given in [10] (Eq. (3.27)) for the square of 
the partition function Zjjasm of U-turn symmetric ASMs of [ '' ]. The latter is known to co- 
incide with XY2rnizi, . . . , Z2m) [-1] and so we reproduce the first factor of the l.h.s. of Eq. (54). 
More precisely, we find A2miwi, . . . , W2m) = 3^""' niri(l + Zi + z:[^)~'^"'+'^XY2mizi, . . . , Z2mf. 
The odd case can be reduced to the even case by sending one of the Zi to zero (something 
which did not make sense in the elliptic setting), so that for both parities we have 

Aniw,, ...,Wn) = 3^L"/2JL(".+i)/2j -Q^i + + z-Y'-^'xyAzi, • • • ,^n)^ 

1=1 

or in terms of the original quantities, A„(xi, . . . ,x„) = 3^^L'^/2JL("-i)/2Jt^^^(^~;_i^^ ^ _ ^ ,-2„)^. 

The second factor is simply obtained by noting that w = J2 = —1/2 corresponds to 

z = u = e*'^/^, so 

n 

Bn{w„ ...,Wn) = 2-32Ln/2+lJL(n+l)/2j "Q^^ ^ ^ ^ _ _ _ ^ 

i=l 

or B„(xi,...,a;„) = 3-2K2JL(n+i)/2j^^^^^(^^^...^^^^^)2_ pinally, 

= n ( Y = 3-"' + Z, + Z.'fxY^iz,, ZnfxY„^, {Zu . . . , Z^, 

i=l ^ ^ i=l 

which is consistent with Eqs. (53) and (54). 
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A. 2. Determinants. Similarly, the determinants simplify as C 0. Noting that w = 
J3 and w = both correspond to 2; = 0, we conclude that there are only two distinct 
determinants for each parity; Tsuchiya's determinant [30, 20] is known to be equal at a 
cubic root of unity to the symplectic character introduced above [ ] 

n l<i<m, {zf + ZiZj + Zj){l + ZiZj + zfZj) 

m+l<j<2m det 



n l<i<j<m, {Zj - Zi){l - ZiZj) l^i<m {zf + Zi.Zj + z]){l + ZiZj + zf z]) 

m+l<i<j<2m m+l<]<2m ■> ■> 

= XY2mi.Zu . . . , Z2ni) 

and then we have: 

2m 

H2m{Wu . . .,W2m) = S'^^""^) + Z, + Zi\'''+\Y,^{zi , • • .,22m) 

i=l 
2m- 1 

H2m{Wl, W2m-l, J2) = S^^""^) \{ {I + Z, + Z-^)-^+\y,,M. • • • , 22m-l, C^) 

i=l 
2m- 1 

H2m{w,, W2m-U h/^) = ^"^^"^''^ 11(1 + ^^ + ^i'T"^^' XY,^., (^l, . . . , 22m-l) 

i=l 

2m 



H2m+2{Wi, . . . , W2m, ^2, J3/4) = T'^^+^h-^ + Z, + Z'^Y^XY^^^r (^1, • • • , 22m, ^) 

i=\ 

2m 

H2m+2{W,, W2m, ^3, ^4) = 3™("*+l) + Z, + Z^r^^'XY^Jz,, Z2m) 

i=l 
2m-l 

H2m+2{Wu . . .,W2m-U ^2, ^3, J4) = 3'"('"+l) JJ (1 + 2:^ + Z'Y'^XY^Jzi, . . .,Z2m-l,Uj) 



i=l 



A. 3. More determinants. The expression (39) of S2m as a Slater determinant reduces to 
the numerator of our definition of the symplectic character XY2m. (since kj — ^m;m+i— j ~l~ 
2m - j + 1, j = 1, . . .,2m) 

S2m(2l, . . . , 22m) = dct {z^^ — Z^ ^) 
t,j=l,--;2m 

The differential equation (49) reduces to 

2m y d \ ^ 

^ ( ^'dz ) • • • ) ^2m) = m(6m^ - l)S2m(2l, • • • , 22^) 

j=l ^ 

A.4. Homogeneous limit. Finally, aI^^ = H2m = 3~'"^™~^^Xy2™(l, • • • , 1) = 1, 1, 3, 26, 646 .. . 
is the number of Vertically Symmetric Alternating Sign Matrices of size 2m + 1 (also, the 
number of Off-diagonally Symmetric Alternating Sign Matrices of size 2m, and the number 
of Cyclically Symmetric Transpose Complement Plane Partitions of size m with a triangular 
hole cut out), while Ali^_^ = i/2m(^3/4) = 3-('"-^)'xy,„_,(l, . . . , 1) = 1, 2, 11, 170 .. . is the 
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number of Cyclically Symmetric Transpose Complement Plane Partitions of size m. Note 
that the square of the number of VSASMs also appears in the observations of [24]. 

The sequence of numbers = 2'^H2m{J2, Js/a) = 3-™("^-i)xy2™+i (1, • • • , 1, = 

1, 5, 66, 2431 . . . appears as one of the factors of the enumeration of UUASMs in [20]. The 
last sequence, 2-(52„_i/3)i/2 = H^^iJ^) = 3-("^-i)\y,^(l, . . . , 1, u;) = 1, 7, 143, 8398, ... is 
the number of ASMs of order 2m + 1 divided by the number of VSASMs of size 2m + 1. 



Appendix B. The C 1 limit 

Besides the C ~^ limit, there is another trigonometric limit, namely C — ^ 1 or a — )• 0. It is 
expected to be somewhat trivial since the correponding Hamiltonian is the Ising Hamiltonian 
with interaction a^a^. Indeed, we find that the building block of the partition function 
becomes: 



l<i<m 

TT \\ _ m+l<j<2m -*- _ om(m-l) 

J^2m[Wi, . . . ,W2m)\(:=l — ] " , Clet ^ — 2 



Yl (1 - AWiWj) 

i<m 
Cj<2m 

l^* ^3) m+l<j<2m ^ 



l<i<j<m 
or 

m+l<i<jr<2m 

This formula is valid as long as the Wi stay finite as C — ^ 1- One special case is if one Wi is 
equal to J4 = 1/(1 — C)- Then we find instead 

SO that 

X„,(. . .) = 2"("+i)+i 

This is compatible with a constant value of '^n,a = 2"^""^)/^ since X2m = ^'^"'^^'^Im a- 



Appendix C. Differential equations 

We provide here analogues of Eqs. (51) and (52) when H2m (that is, the function H2m 
with all arguments set to zero) is replaced with H2m{S), S C {J2, J3, J4} (again, with all 
other arguments set to zero). Because of the permutation symmetry w.r.t. {J2, J3, J4}, we 
only need to provide one formula for each possible cardinality of S. When taking derivatives 
w.r.t. u or V, the convention is that the arguments that are specialized to J2, J3, are among 
the u's. 

After transposition (for display purposes), Eq. (51) is of the form 

/ 92 Q2 Q2 Q2 d d 

\ouov ow^ oa'^ ouoa ou oa 
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For itself, the matrix P is 








2(4m+l) 






-2m(a-l)(a+8) 
-2(-4m-l) 




-2Tn{a~l){a+S) m^Qm^+2m^— am— 2m— 40— 14^ 



2{2m+l) 
4m 


~2m{a~l){a+8) 
m(m^ {a+2) 




2(4m2+m-2)a 
— 2m(4m+l)a 



m(2m+l)(a-4)2 




\ (m-l)m2(a+2) 



(m— l)m 



-(m— l)m^ (a— 4) — (m— l)m^ (2m+l)(o— 4)o / 



For H2M2] 



( 



2(2m+l) 

4m 

-2(m-l)(a-l)(«+8) 

-2(l-4m) -2m(a-l)(a+8) 

2(4m— 1) m(^am?+2m'^—ain—4in+a+4^ m{2mo?—o? — l&ma+Z2m}j 

-2(m-l)(a-l)(a+8) (m-l)(am2+2m2-am-4m-4«-12) 

(m— l)^m 



2(4m2-5m-l)a \ 
— 2m(4m,— l)a 





\ (m-l)2(am+2m-2) 



'{m—l)m(am—4m+A) — (m— l)m^a(2am— 8m— «+8) / 



For H2miJ3, Ji)- 



2(2m+l) 

4m 

-2(m-2)(Q-l)(a+8) 

2(4m-3) -2m(o-l)(«+8) 

2(4m-3) ni{^am? +2m? -am+a+2) m{2ma'^ -a^ -l&ma+l&a+32m-i2^ 

-2(m-2)(a-l)(a+8) {;m-2){^an? +2m? -am-ia-l-i) 

(m— 2)(m— l)m 



2(4m2-llm+4)a \ 
-2m(4m-3)a 





\ (m— 2)m(am+2m— a) 



-(m—2)m(am—Am—a) —{m—2){m—l)ma{2am—8m—a) / 



For H2m{J2, J3, Ja)- 



/ 2(2m+l) 

4m 

-2(m-3)(a-l)(a+8) 

2(4m-5) -2m(«-l)(a+8) 

2(4m-5) m(m2-m+l)(a+2) 

-2(m-3)(a-l)(Q+8) {m-3)(^am'^+2m'^-am-2m-4a-14^ 



2 m— l)(4m-13)a \ 
— 2m(4m— 5)q! 



m(2m-3)(a-4)2 




\ (m-3)(m-l)m(a+2) 



(m— 3)(m— l)m 



-(m-3)m2(a-4) -(m-3)m2(2m-3)(«-4)a / 



As to Eq. (52): 



Co H2( m+l) (m— 1 ) — Ci H2mH2m ~ C2 (-f^2m) + ^3 -ff2m-f^2m + ^4 



2 

2m 
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The coefficients for i72m are: 

Co = 4Q;(4m - l)(4m + l)2(4m + 3) 

Ci = 2{a - ifaia + 8)^(4m + if 

C2 = 4(q; - l)^a{a + 8)^(4m - l)(4m + 3) 

C3 = (a - 1)(q; + 8)(q;^ + 28q; + 2Wm^ + 304Q;m^ - 256m^ + 20a'^m + 20Sam - 192m - 32) 
C4 = 60;^ - 24q; + 4Q;3m^ - 1008Q;^m^ + 3408Q;m^ + 512m^ - 4Q;^m^ - 984Q;^m^ + 4032Q;m^ 
-128m^ + - \^2o?n? + 1028Q;m^ - 320m^ - Q?m + 28Q;^m - 44Q;m - 64m 

For E^mW: 

Co = 4Q;(4m - 3) (4m - l)^(4m + 1) 
Ci = 2(4m - + afa{% + af 

C2 - 4(4m - 3)(1 + 4m) (-1 + afa{^ + q)^ 

C3 = (-1 + a)(8 + a) (64m - 256m^ - 4q; - 96mQ; + 304m^Q; + - \mo? + 2Am'^a^) 
C4 = -128m + 768m2 - 1152m3 + 512m^ - 36a + 336ma - 204m2a - 2784m-^a + 3408m^a 
+18a^ - 126ma^ - 6m^a^ + 1032m^a2 - 1008m^a^ - ma^ + Qm^a''^ - 12m^a''^ + Am^a^ 

For //2m (^3, ^4): 

Co = 4a(4m - 5) (4m - 3)2(4m - 1) 

Ci = 2(4m - 3)'\-l + a)^«(8 + 

C2 = 4(4m - 5)(4m - 1)(-1 + afa{8 + 

C3 = (-1 + a)(8 + a)(-192 + 448m - 256m^ + 204a - 512ma + 304m^a + 2W 

-44ma^ + 24m^a^) 

C4 = 384m2 - 896m^ + 512m^ + 720a - 5208ma + 11892m2a - 10848m^a + 3408m^a 
-90a^ + mima^ - 2958m'^a^ + 3000m^a2 - 1008m^a^ + 3ma^ - 3m2a^ - 4m^a^ + 4m^a^ 

For H2m{J2: Js: Ji)- 

Co = 4a (4m - 7) (4m - 5)2(4m - 3) 
Ci = 2(4m - 5)^(-l + a)2a(8 + af 
C2 = -4(4m - 7) (4m - 3)(-l + a)2a(8 + af 

C3 = (-1 + a)(8 + a) (-480 + 704m - 256m^ + 540a - 816ma + 304m^a 

+45a^ - 68ma^ + 24m^a^) 
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C4 = -960m + 2368m2 - 1920m^ + 512m^ + 8400a - 27740ma + 33572m^a 

-17664m^a + 3408m^a - 2100a^ + 7240ma^ - 9142m^a^ + 5016m^a^ 
— 1008m'^a;^ — 5ma^ + 13m?a^ — 12m^a^ + Am'^a^ 
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